% 创建时间：2025-06-20 23:29:12
% 作者：兪馳
clc;
clear;

A=[8,1,6;
    3,5,7;
    4,9,2]; 

Q=[16,4,1;
    9,3,1;
    4,2,1];

X=dlyap(A,Q), 
norm(A*X*A'-X+Q)  % 精度验证